Anomalous Josephson current through a ferromagnetic trilayer junction 
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We studied the anomalous Josephson current appearing at zero phase difference in junctions 
coupled with a ferromagnetic trilayer which has noncoplanar magnetizations. A tt/2 junction with 
an equilibrium phase difference ir/2 is obtained under suitable conditions. The equilibrium phase 
difference and the amplitude of the supercurrent are all tunable by the structure parameters. In 
addition to calculating the anomalous current using the Bogoliubov-de Gennes equation, we also 
developed a clear physical picture explaining the anomalous Josephson effect in the structure. We 
show that the triplet proximity correlation and the phase shift in the anomalous current-phase 
relation all stem from the spin precession in the first and third ferromagnet layers. 

PACS numbers: 74.50.+r, 74.45.+C, 74.78.Na 
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INTRODUCTION 



Usually the supercurrent in a Josephson junction van- 
ishes, when the phase difference between the two super- 
conductors is zero, and in the tunneling limit the current- 
phase relation (CPR) is sinusoidal I(ip) = I c sm((p)J* 
Recently some studies 2 - - — found an anomalous Joseph- 
son current flow I a exists even at zero phase difference 
(tp = 0). The anomalous supercurrent is equivalent to 
the presence of an additional phase shift ipo in the con- 
ventional CPR, i.e., I((p) — I c sin(ip + ipo). In fact, such 
CPRs have been predicted for Josephson junctions of un- 
conventional superconductors^ - — but the experimental 
verification is still lacking. Recent studies have shown 
that the anomalous supercurrent can also exist in junc- 
tions with conventional s-wave BCS superconductors if 
both spin-orbit interaction (SOI) and a suitably oriented 
Zeeman field are present in the coupling layer.™ - -- These 
studies revealed that the anomalous effect in conventional 
junctions has some intricate physics. More interesting, an 
anomalous Josephson current can also appear in super- 
conductor (S)-ferromagnet(F) hybrid structure without 
SOL— £ In Grein's study, 5 a SFS hybrid structure with 
two spin-active interfaces was considered. The two spin- 
active interfaces are critical to the triplet proximity effect 
and the anomalous supercurrent in the structure, but the 
physics is still unclear. 

In this study, we generalize the two spin-active inter- 
faces to two ferromagnetic layers with finite thicknesses 
and clarify the physical mechanisms responsible for the 
anomalous supercurrent. In such SFFFS structures, we 
find that the triplet proximity correlation and the phase 
shift in the anomalous CPR all stem from the spin pre- 
cession in the first and third F layers. According to the 
symmetry analysis^ an anomalous supercurrent is pos- 
sible when the symmetries of the time-reversal opera- 
tor T and its combination with a spin rotation operator 
with respect to an arbitrary spin quantum axis n er n T 
are broken at the same time. As a result, the simplest 
superconductor (S)-ferromagnet (F)-superconductor (S) 
junction for achieving an anomalous Josephson current 



requires the F layer to be a ferromagnetic trilayer with 
noncoplanar magnetizations for breaking the symmetry 
of the operator a n T. SFFFS junctions where the mag- 
netizations of the three ferromagnetic layers need not be 
noncoplanarii^ and SFS junctions with inhomogeneous 
magnetization^ - — , have been studied in order to under- 
stand the effects of triplet correlation induced in the F 
layers. Controllable 0-tt transition and spin-triplet super- 
currents have been realized experimentally recentl y 21 ' 22 . 
In our study, we found that triplet correlation is also an 
important condition for the anomalous supercurrent 32r— 

We consider a junction consisting of two conventional 
s-wave superconductors coupled by a ferromagnetic tri- 
layer with noncoplanar magnetizations. For convenience, 
hereafter we denote the three F layers sequentially by 
Fi,F 2 ,F 3 . We start with the typical situation where the 
magnetizations are along the x, y, z axes respectively (i.e. 
an SF x F y F z S junction), as shown in the upper panel of 
Fig. [T] This junction is a n/2 junction with an equi- 
librium phase difference n/2 under suitable conditions. 
The equilibrium phase difference can be tuned by the 
lengths, the exchange energies, and the magnetization 
orientations of the F± and F3 layers. And the amplitude 
of the supercurrent can be tuned by the barriers between 
the F layers or by the length and the exchange energy 
of the middle F2 layer. In this regime the Josephson 
junction can also act as a supercurrent rectifier i 28 i 29 

The paper is organized as follows. In Sec. II we present 
the model and solve the scattering problem for quasi- 
particles based on the Bogoliubov-de Gennes equation. 
The Josephson current and Andreev bound states can be 
obtained from the scattering matrices. In Sec. Ill we 
show the numerical results for the anomalous supercur- 
rent and corresponding Andreev bound states and reveal 
the physics. A conclusion and remarks will be given in 
Sec. IV. 



II. MODEL AND FORMALISM 

In the numerical calculation, we consider SF1F2F3S 
junctions with various lengths and exchange energies for 
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FIG. 1: (Color online) Upper panel: Schematic diagram of 
the SF x F y F z S junction where two barriers are present be- 
tween the F layers. Lower panel: Schematic illustration of 
the formation of Andreev bound states with triplet correla- 
tion in the F y layer due to spin precession of electrons and 
holes in the F x and F z layers. 



each F layer and various barrier strengths for the two 
barriers between the F layers. The transport direction 
is along the x axis. The three F layers have the thick- 
nesses, Li,L2,Li3, the exchange energies, ht,h2,ha, and 
the magnetization orientations, </>i), (#2, <fo), (#3, ^3) 
in spherical coordinates. The effective Hamiltonian of 
the system is given b y 30 ' 31 
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where e k = ^k(^x + ^y ~ k F ) + ^ with ^f the Fermi 
wave number, U = Uq[S(x — L{) + S(x — L\ — L 2 )] 
represents the two barriers between the F layers, and 
h z = hcosO, h xy — hsmde 1 ^, with h the strength 
and (8, <f>) the orientation of the exchange field; A(x) = 
A[Q{-x)e lip / 2 + <d{x-L)e- lip /' 2 ] describes the pair poten- 
tial with L = L \ + L2 + L3 and A the bulk superconduct- 
ing gap and tp = ipi — ipx the macroscopic phase differ- 
ence of the two superconductor leads. The temperature 
dependence of the magnitude of A is given by A(T) = 
A(0) tanh( 1. 74 t/Tc/T- 1)32 

Since the transversal mo- 
mentum components are conserved and not important 
to the total Josephson current, we consider the ques- 
tion in the one-dimensional regime for simplicity. The 
Bogoliubov-de Gennes equation can be easily solved for 
each superconductor lead and each F layer respectively. 
The scattering problem can be solved by considering the 
boundary conditions at the interfaces. Each interface 
gives a scattering matrix. The total scattering matrix 
of the system can be obtained by the combination of all 
these scattering matrices of interfaces. From the total 
scattering matrix, we can obtain the Andreev reflection 



amplitudes ai CT and a^a of the junction where a\ a is for 
the reflection from an electron-like to a hole-like quasipar- 
ticle and a^a is for the reverse process with <j representing 
the spin. The stationary Josephson current can be ex- 
pressed in terms of the Andreev reflection amplitudes by 
using the temperature Green function formalism^ 3 . 
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where k n , a\ an , and 02 



are obtained from fc+, fc~ , 
ai a , and a^a by analytic continuation E — > iuj n . kf is the 
wave vector for electron or hole in the superconductors 
and the Matsubara frequencies are w„ = 7rfcsT(2n + 1), 
n = 0, ±1, ±2, • • • , and fi n = + A 2 . 

The discrete spectrum of the Andreev bound states can 
be determined by using the condition 34 



det[I - R 2 PR 1 P] = 



(3) 



where Ri,R2,P are 4x4 matrices, P is the propaga- 
tion matrix of modes in the F2 layer, and R\ (R2) is the 
reflection matrix of the right-going (left-going) incident 
waves. 

In order to study the spin properties of the Andreev 
bound states formed at F2 layer, we can also work out the 
Green's function G(x, x' , E) in F2 layer which is a 4 x 4 
matrix. 3 - 5 Now it is convenient to take the eigen spinors 
of F% layer, i.e., spin-parallel and spin-antiparallel with 
respect to the exchange field h 2 as the unit vectors of 
the spin space. Then the spin current in F% layer can be 
evaluated by^S 
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where 7_|_ (7_) is the charge currents of electrons with 
parallel spin (antiparallel spin) and obviously satisfies 
h = /++/-. 



III. RESULTS AND DISCUSSION 

We start with the typical noncoplanar magnetization 
configuration, i.e., the SF x F y F z S junction. For simplic- 
ity, we introduce the dimensionless units: the energy 
E EEf, the wave vector k —> kfc^, the coordinate 
x — > x/fcp, and the strength of exchange field h — > hEp. 
All physical quantities are expressed in the dimension- 
less units in the rest of the paper. The superconductors 
considered are characterized with A = 10~ 3 which corre- 
sponds to the BCS coherence length at zero temperature 
£0 = 2/ttA k, 636.6. 
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FIG. 2: (Color online) Charge and spin currents 7+, 7_, I e , 
I 3 versus p for the SF x F y F z S junction. The barrier strength 
Uo — for (a) and Uo = 2 for (b). The strength of exchange 
fields hi = h% = hz = 0.05, and the lengths of F layers 
Li = L 2 = L 3 = 10tt. The temperature T/T c = 0.5 with T c 
being the critical temperature. 

Fig. [2] shows the charge and spin currents 7+, J_, I e , I s 
as functions of the phase difference cp for the SF x F y F z S 
junction. The corresponding Andreev bound states are 
shown in Fig. [3l It is interesting to note that when there 
is a barrier between the F layers, (Uo ~ 2), there is a 
significant anomalous Josephson current. When there is 
no barrier Uq = 0, the anomalous Josephson current is 
nearly zero. This interesting dependence on the barrier 
strength Uq will be explained below in terms of the spin 
characteristics of the Andreev bound states in the F y 
layer. 

Firstly, it is useful to point out a large spin current 
exists in the F y layer, implying that the superconduc- 
tivity correlation is mainly triplet in the F y layer. This 
is easily understood by considering the formation of an 
Andreev bound state in the F y layer. A right-going elec- 
tron with spin parallel to the y axis (1,«) T from the F y 
layer will have its spin precessing about the z axis in the 
F z layer before it reaches the right superconductor. Af- 
ter the Andreev reflection from the right superconductor, 
a hole with reverse spin goes left and its spin continues 
to precess. The one-way angle of precession is approxi- 
mately (fc+ — fc_)L3 h^L^ where fc+ (fc-) is the wave- 
vector of up-spin (down-spin ) quasi-particle. Thus if 
the condition /13L3 = nix + n/2 (n is an integer) is sat- 
isfied, the reflected hole from the right superconductor 
will have its spin parallel to the incident electron's spin 
in the F y layer. An Andreev bound state is formed, after 
this reflected hole travels through the F y and F x layers 
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FIG. 3: The energy levels of the Andreev bound states Eb- 
The strength of exchange fields hi = h% = hz = for (a) and 
hi = h% = /13 = 0.05 for (b), (c), (d). The barrier strength 
Uo = for (b), U = 1 for (c), and U = 2 for (d). The other 
parameters are the same as those in Fig. [2] 



and Andreev reflected from the left superconductor and 
changes to an electron to move right to finish a cycle. If 
the spin rotation angle in the F x layer satisfies the same 
condition h\L\ = mr + 7r/2. The electrons and holes 
have identical spins (parallel to the y axis) in the F y 
layer and the Andreev bound state formed has complete 
triplet correlation in the F y layer, as schematically shown 
in the lower panel of Fig. [T] Triplet correlation can exist 
in other different type of magnetic inhomogeneity too. 
Bergeret et al^Z- have studied S/F/S junctions with spi- 
ral magnetization in the F layer and found spin triplet 
correlation there. In the present model, we found two 
Andreev bound states below the Fermi level with com- 
plete triplet correlation; one is "spin- up" (with respect 
to the y axis), which carries the current J +) and one is 
" spin-down" , which carries the current i_ , as shown in 
Fig. [3J In the short junction limit, the Josephson current 
is totally carried by the Andreev bound states^ 

Besides complete triplet correlation in the F y layer, 
another interesting feature noted in Fig. [5] is that J+ has 
a phase shift of 7r/2 while 7_ has a phase shift of — ir/2 
compared with the conventional CPR. So, these two 
currents move in opposite directions. Now we follow the 
Andreev reflection processes occurring in the formation 
of the bound states to find out the phase shift. For 
simplicity, we assume hi = h% = ha = h ; thus, the 
wave-vectors of "spin-up" (+) and "spin-down" (— ) 
electrons (holes) with energy E at each F layer are 

k± h) = \Jk 2 F + p e (h)E T h with p e[h) = +(-)l. In the 
short junction limit and the limit E <C h <C Ep , we have 
k± s» k± as k± — kp =F We start with a right- going 
"spin- up" electron at the position x = Li + 0, the wave 
function can be written as (l,i,0,0) T . The electron 
moves right and acquires a phase e zk+L2 when it arrives 
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at the interface x = L\ + L2. To simplify the discussion 
we focus on the Andreev reflections at the F/S interfaces 
and ignore the normal reflections at the barriers which 
affect only the amplitude of the supercurrent but not 
the phase shift. When the electron travels through 
the F z layer, its spin precesses. The state becomes 
(e ifc + i3 ,ie lfe - i3 ,0, 0) T e ik + L2 when the electron arrives 
at the interface x = L\ + L2 + £3. Then, the electron 
is reflected as a hole with reverse spin and the hole 
wave function is (0, 0, -ie lk - Ls , e ik+L 3 y e ik+L 2 ^_ e i<p/2 

where u = + fl/E)/2, v = ^(1 - tt/E)/2 with 
il = y/E 2 — A 2 . The algebraic derivation is not 
shown here for space limitation and the approximation 
kj_ ps k\ ps k%_ ps kp has been used in the deriva- 
tion where k s + (fcl) is the wave- vector of electronlikc 
(holelike) quasiparticle in the superconductors. The 
Andreev-reflected hole moves left and has its spin 
rotated in the F z layer again and then goes back 
to the F y layer x = L\ + L2 — 0. Now the wave 
function becomes (0, 0, -ie lhL \ e - ihL *) T e ik + L * % e lv ' 2 = 
(0, 0, 1, -i)T e ik+L 2 ji e iip/2 where the con dition hL 3 = it/2 
has been used. The wave function describes a "spin- 
up" hole with respect to the y direction. Then the 
hole goes left through the F y layer and acquires 
a phase e~ lfe + 2 . So the wave function becomes 
(0, 0, 1, — i) T ^e lip / 2 when the hole arrives at the interface 
x = L\. Consequently, the hole has its spin precessed 
in the F x layer and moves left to the interface x = 
with the wave function ± [(1 - i)(0, 0, 1, l) T e _lfc + Ll + 
(1 + i)(0,0,l,-l) T e~*- Ll ] ie 1 ^ 2 . The hole is 
Andreev-reflected as an electron with reverse 
spin described by \ [(1 - i)(l, -1, 0, 0) T e~ lk + Ll - 

(l + i)(l,l,0,0) T e- ife - il ] {^fe lv . Then the electron 
goes through the F x layer again and back to the starting 
position x = L± + to finish a cycle. The final wave 
function is (1, i, 0, 0) T e l7r / 2 e lv where hLi = tt/2 
is used. Comparing with the initial wave function 
(l,i,0,0) T , we can see the phase shift of the "spin-up" 
Andreev bound state is indeed tt/2 when considering a 
conventional CPR. In the same way, we can find out the 
phase shift of the "spin-down" Andreev bound state is 
—tt/2. In this round-trip cycle of the quasi-particle, we 
can clearly see that the phase shifts in Andreev bound 
states come from the spin precession of electron and hole 
in the F x and F z layers. 

If we neglect the second-harmonic term in the CPR, 
the charge current carried by the two Andreev bound 
states can be written a.sM- 

7+«7°shn>+|), 7_«7°shn>-|), (5) 

where 7? (7°_) is the amplitude of the "spin- up" ("spin- 
down") charge current. When the barriers are absent, 
the normal scattering at the two F/F interfaces can be 
ignored and we can have 7? ss 7°_. As a result, the total 
charge current I e = I + + 7_ is very small and only the 
second- harmonic term remains, as shown in Fig. [5] (a). 
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FIG. 4: (Color online) J e versus <p for the SF x F y F z S junction 
with different 83: varying 63 from to n with a step 7r/4, 
4>3 = 0. The other parameters are chosen as: Uo = 2, hi = 
h 2 = h 3 = 0.1, Li = L 2 = L 3 = 5tt, T/T c = 0.5. 



At zero phase difference, the charge current is very small 
and the spin current in the F y layer is almost a pure spin 
current. 

When the barriers are present, the normal scattering at 
the barriers reduces the amplitudes of the two charge cur- 
rents and J° • The transmission probability through 
the double delta function barriers of electrons or holes de- 
pends on the wave-vector of the particle in the F y layer 
and reaches the maximum when resonance transmission 
occurs. Here in the F y layer, the "spin-up" Andreev 
bound state couples a "spin- up" electron with a "spin- 
up" hole which have the same wave- vector k + ss kp — k 
while the "spin-down" Andreev bound state couples a 
" spin-down" electron with a " spin-down" hole which have 
the same wave- vector fc_ ss kp + ^. The difference in 
the wave vector between the two Andreev bound states 
leads to the difference in the transmissions through the 
F y layer. Consequently, we can make a large difference 
between 1+ and 7°_ as shown in Fig. [2] (b) by using two 
barriers as well as suitable exchange field strength and 
length of the F y layer. In this way, an anomalous Joseph- 
son current appears at zero phase difference. The CPR 
of the junction has a phase shift of ±ir/2 in comparison 
with the conventional CPR where the sign of the phase 
shift depends on the relative magnitude of 7+ and 7°_ . 

Since the phase shift of the anomalous CPR stems from 
the spin precession of electrons and holes in the F x and 
F z layers, we can modulate the phase shift by tuning 
the parameters of these two layers. If the conditions 
0i = 4>3 — and h\L\ = I13L3 = (n + l/2)7r are satisfied, 
the complete equal-spin triplet correlation in the F y layer 
is maintained. Now the phase shifts of the "spin- up" and 
"spin-down" Andreev bound states are + O3 — 0\) 
according to the above discussion. Fig. 2] shows the 
tuning of the equilibrium phase difference by varying 63. 
Bergeret et al~± have found that the relative orientation 
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FIG. 5: (Color online) I e versus ip for the SF x F y F z S junction 
with different Uo and /12: (a) varying Uo from —2 to 2 with 
a step 1, hi = 0.1, Z/2 = 5ir; (b) varying ^2 from 0.05 to 0.15 
with a step 0.02, Li = 10-7T, Uo = 2. The other parameters 
are chosen as: hi = /13 = 0.1, Li = L3 = 57r, T/T c = 0.5. 



FIG. 6: The anomalous supercurrent at zero phase difference 
I e (<p = 0) for the SF x F y F z S junction as functions of: (a) 
hi with 1/2 = 10-7T, and (b) Li with hi = 0.1. The other 
parameters are chosen as: hi = /13 = 0.1, L\ = L3 = 57r, 
C/ = 2, T/T c = 0.5. 



of the two magnetizations in S/F/I/F/S junctions can 
change the critical current. For the particular structure 
considered in Fig. 01 which is different from theirs, the 
orientation of the magnetization in the third layer has no 
strong effect on the critical current. However, for other 
values of /i 2 , the orientation can also modify the criti- 
cal current. On the other hand, the amplitude or even 
the sign of the supercurrent can be changed by the bar- 
rier strength, the exchange field strength or the length 
of the F y layer, as shown in Fig. [5] The dependence 
of the supercurrent on the barrier strength is because of 
the condition of resonance transmission through double 
delta barriers sin(2kLi) = —4U / (£/ 2 +4) with k the wave 
vector of particles. Fig. [6] shows the anomalous super- 
current at zero phase difference as functions of hi and 
Li for the SF x F y F z S junction. It is noted that the de- 
pendence on hi exhibits a period of 2ir/Li = 0.2 which 
confirms the occurrence of resonance transmission of elec- 
trons and holes through the F y layer. And the depen- 
dence on Li exhibits two periodic behavior. One period 
is nearly ir and the other is 107r. Because the wavevec- 
tor of "spin-up" electrons and holes is k + ss kp — 4^- 
while the wavevector of "spin-down" electrons and holes 
is fc_ kp + 4p in the F y layer, the period n stems from 
27r/2fcp = 7r with kp = 1. The other period 107r stems 
from 2ir/2hi = lOir which is related to the difference of 
wavevectors k + — fc_ = — hi. 

If anyone of the two conditions </>i = </>3 = and 
h\Li — h s L 3 = (n + l/2)n is not satisfied, the pure 
equal-spin triplet correlation in the F y layer is changed. 
For example, if we vary the length or the exchange field 
strength of the F x and F z layers, the spin precession angle 
of electron and hole in a round trip in the F x and F z lay- 
ers is not 7T any more. The reflected hole will have both 
the same spin component and the opposite spin compo- 
nent to the incident electron in the F y layer. Now the 
correlation is the mixing of singlet and triplet. But only 
the triplet correlation can contribute to the anomalous 
Josephson current, so the anomalous supercurrent is re- 
duced with increasing singlet component. Fig. [7] shows 
that both the amplitude of the supercurrent and the equi- 
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FIG. 7: (Color online) I e versus tp for the SF x F y F z S junction 
with different Li: varying L\ from to 57r with a step 7r, 
Lz — Li. The other parameters are chosen as: Uo = 2, 
hi=h 2 = h 3 = 0.1, Li = 5tt, T/T c = 0.5. 

librium phase difference is tuned by the length of the F x 
and F z layers. To study the characteristics of Cooper 
pairs in the F y layer in detail, the pair function can be 
defined by the anomalous Green function and be decom- 
posed into four component s 23 ' 25 

3 

^2G*{x,x)=i^2Mx)a v oi, (6) 

w„>0 v=0 

where G" 1 is the anomalous electron-hole correlation 
function, cto is the unit matrix and a v (y = 1,2,3) are 
three Pauli matrices. In Eq. ([5]) , the frequency summa- 
tion is only made over positive frequencies because the 
triplet pair functions are odd functions of frequency, fo 
(/a) is the pairing function of spin-singlet (spin-triplet) 
pairs with spin structure of - (+) lit)] /V%- The 
pairing functions of |tt) an( l pairs are given by 
/ft = if 2 - fi and f u = if 2 + f u respectively. Fig. 
[5] (a) shows the absolute values of pairing functions at 
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FIG. 8: (Color online) The absolute values of pair functions 
in the F y layer (a) and the anomalous supercurrent at zero 
phase difference (b) as functions of the length of the F x and 
F z layers for the SF x F y F z S junction. Lz = L\. The other 
parameters are the same as in Fig. [7] 



the center of the F y layer as functions of the length of 
the F x and F z layers for the SF x F y F z S junction. The 
equal-spin pair functions and opposite-spin pair functions 
oscillate with the length of the F x and F z layers which 
determines the angle of spin precession of quasiparticles 
in these two layers. Compared with the anomalous super- 
current shown in Fig. [5](b), we can see that the anoma- 
lous supercurrent is nearly proportional to the equal-spin 



triplet correlations. 



IV. CONCLUSION 

In summary, we predict a tunable anomalous Joseph- 
son effect in SF1F2F3S junction where the three F layers 
have noncoplanar magnetizations. The superconducting 
correlation can be completely triplet in the F2 layer due 
to the spin precession of electrons and holes in the F± and 
F3 layers. If the condition h\Li = /13L3 = (n + l/2)w 
is satisfied, an electron incident to the left (right) su- 
perconductor will precess its spin by ^ in the F\ (-F3) 
layer before it arrives at the supercondcutor and the 
Andreev-reflected hole proceeds to precess the spin by 
^ when it goes back to the F2 layer. Thus the Andreev- 
reflected hole will have the same spin with the incident 
electron and the complete triplet correlation arises in the 
F 2 layer. The two spin-resolved Andreev bound states 
carry two spin-polarized supercurrents which have oppo- 
site phase shifts and different amplitude thus leading to 
an anomalous Josephson current. And the phase shift in 
the anomalous current-phase relation is also a result of 
the spin precession of electron and hole in the F\ and F3 
layers. The equilibrium phase difference of the anomalous 
supercurrent can be tuned by the lengths, the exchange 
energies, and the magnetization orientations of the F\ 
and F3 layers. And the amplitude of the supercurrent 
can be tuned by the barriers between the F layers or by 
the length and the exchange energy of the F2 layer. 
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